A method to determine the polarization of high energy gamma - rays 
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A method to determine the polarization direction of a very high energy (50 MeV - 30 GeV), 
linearly polarized gamma ray is analyzed and discussed. We study the differential cross section for 
pair production and show that the only relevant piece of information comes from the reconstruction of 
the pair trajectory. In particular the polarization information is coded in the azimuthal distribution 
of the created pair. We show that the Azimuthal Ratio (AR), defined as the ratio of produced 
pairs parallel to the polarization vector to pairs perpendicular to this vector, yields a very practical 
method to determine incoming polarization. We investigate the behavior of this function for different 
' angular resolutions of the measuring apparatus and different energies of the incoming pairs. 
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I. INTRODUCTION 



PsJ ' The study of high energy gamma rays produced by astrophysical sources give us a unique window to study very 
C^l unusual objects in our universe. The GRO [Q has provided plenty of new data to study gamma ray bursts, and to 
obtain detailed maps of our galaxy in several energy ranges At present there are several projects to improve 
. the resolution of GRO |^J^ and give us a more accurate description that will enable us to decide if the gamma ray 
■ sources are extended or compact objects. It is worth mentioning that neither GRO nor the new projects have been 
' designed to detect the polarization of those rays. However, one knows that polarized gamma rays are very important 
, in astrophysics since they constitute a distinctive byproduct of matter accretion by a black hole or neutron star. 
\C ' Thus, it could be desirable to measure incoming polarization together with the other variables and perform a more 
comprehensive study of those astrophysical processes. 

We address here the problem of measuring the polarization of gamma rays in the energy range 50 MeV - 30 GeV. 
\ The interaction of gamma rays with matter in this range is dominated almost exclusively by pair production. Thus, 

• the detection of those rays is based on measuring the electron - positron pair created by these photons (Fig.|l|). 

Oh' In this work we extensively study the cross section of pair production looking for what characteristics of this process 
Q , can be useful to determine the polarization of the gamma - ray. Since most astrophysical processes of interest generate 
$H ■ linearly polarized rays we well assume that the incoming ray is polarized along a definite direction (along the x - 
c/3 direction in Fig.|^). 

• We first show that it is not possible to determine the incoming polarization by measuring the spin of the resulting 
^ ] pair. This leaves the study of the particle track as the only mean to detect the polarization direction. In particular, 

the polarization information lies in the azimuthal distribution of the particles moments so we calculate this probability 
distribution integrating the differential cross section in energy and polar angles. 
^ I The idea to detect polarization of a gamma ray by measuring the azimuthal distribution of the created pair was 
■ ■ ' brought up in the 50's. It was then determined that the asymmetric ratio (AR) of the number of pairs produced 
in a plane parallel to the polarization vector, to the number perpendicular to that direction is a good parameter to 
describe this distribution |^,^. Although at that time only cross sections for coplanar events were used, Maxinion 
and Olsen later showed that this ratio is also useful when the events make a small but finite angle Ae around the 
coplanar direction [Q. 

However, since the differential cross section for this process is a very complicated function of the relevant angles, 
the above mentioned authors used approximations with different, and sometimes opposite, results. As an example, 
Berlin et al. Q calculated a higher probability of pair production perpendicular to the polarization vector, whereas 
Wick Q obtained the opposite result. Maximon and Olsen ||^, using quantum mechanical calculations, found that the 
azimuthal relation asymmetry depended on the angular resolution with which one measures these probabilities. Thus, 
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FIG. 1. Angles occurring in the pair creation 



depending on this resolution, one might have a higher probability of pair production in either direction. However, 
since their calculations were only done for small deviations from coplanarity it was not clear what would be the 
behavior of AR as a function of the azimuthal angle. 

It is worth mentioning that none of these authors computed the azimuthal probability distribution which is necessary 
to perform Monte Carlo simulations test to the prototypes. 

In this work we calculate this spatial azimuthal distribution for high energy gamma rays (50 MeV - 30 GeV) directly 
from the differential cross section by integrating over the polar angles and the positron (electron) energy. 

The integration is made done using a Monte Carlo method and the VEGAS routine of Numerical Recipes ||^ . To 
check the method, we calculate the total cross section and find that it yields the tabulated values. 

In Section II we first compute the scattering cross section for pair production without summing over spins or 
polarization directions. This would be a standard problem for a course on QED but it is not discussed on any 
textbook on the subject. The main result obtained is that, in the high energy limit, the spinorial part factorizes out 
in the final expression thus showing that the polarization angle does not mix with the spin directions. We then sum 
over the spin directions obtaining the standard May formula 

In Section III we compute the spatial azimuthal distribution by integrating the cross section over the energy and 
polar angles. We also parametrize this distribution with a 'fitting function' of the two azimuthal angles. 

Next, in Section IV, we compute the asymmetric ratio AR as a function of the indeterminacy in the azimuthal 
direction. This is done by integrating over the angle Ae and an appropriate energy range of the produced pair. 
Calculations are repeated for different incoming energies of the gamma ray. 

Finally, in Section IV we argue that the best method of detecting incoming polarization is to select the angular 
resolution that yields the maximum value of AR. 



II. RELATIVISTIC CROSS SECTION FOR LINEARLY POLARIZED GAMMA RAY 



The pair production by a high energy incident photon can be analyzed with a very good approximation by assuming 
the incident photon interacts with an external potential, which is assumed to be a spherically symmetric, static one 
(Born approximation). In the following, Ur and Vs (r, s — 1, 2) denote spinors for the electron and positron respectively, 
such that IFur' = Srr'Ur, TVvr' = (1 — Srr')vr, whcrc 11"'^ (11^) is the positive (negative) helicity projection operator; 
and A+ (A^) is the positive (negative) energy projection operator. That is, 

n'-(p) = i[I + (-l)Vp] , A±(p) = ^(±^ + m) . 
The Feynman amplitude for a pair defined by these spinors is = Aia + ■M-b, with 

and 
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m 



.M2 



{i — a, b) 



The absolute value of Ai, 
can be computed with the aid of the above mentioned projectors. Defining 



Mf = \Maf + \Mb\^ + 2?fi{MaMb) 



(i = a, b) 



|gr[(K)'-m2] 

using the identity Ur{p)ur{p) ~ WA'^{p) and similarly for v, we find that 

\Ma\' = -cl tr[ A+(pi) Ta A-(p2) ri nn , 
?fi{MaMb) = -caCb K{tr[A+(pi)rfcff A-(p2)rinn} . 

Since we are interested in very energetic gamma rays, we can take the ultrarelativistic limit for the projection 
operators, namely, 



(1) 
(2) 
(3) 



2m 2m 



Inserting the above operators in Eqs.(|l|,^,^ and recalling that no linear terms in must appear in the final expressions 
since we are dealing with real quantities, we obtain 

^2 



\Ma\ 



\Mb\'' 



^(MaMb 



16m^ 

CaCb 



Note that the spinor contribution factorizes out in the three equations, and that only two final states (which are 
equally probable) are allowed, independently of photon's polarization. Thus, it is not possible to know the latter from 
the spins of the pair. We now proceed by adding over the two possible final states. Assuming that the photon is 
lineally polarized we obtain 



\M\' = -'^ 



2m?\q\'^ 



{^Pif 



[pikf 



[AEl - |,f ) 



2(pie)(p2e) 



, (P2e) , (pie) 



(4) 



{pik){p2k) (pik) {p2k) {pik){p2k) 

The cross section for pair production by lineally polarized gamma rays in the high energy or ultrarelativistic limit 
is thus obtained by inserting the Feynman amplitude (0) in 



d(7 = (27r)^ -r-^o 
47r 



E+E-Uj\^'^{2my 



\M\'d''pid''p25 {E+ +E^-Lu) 



Computing the scalar products and using natural units with h/2TT — c ^ 1 this gives 
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with 



-2 [E{uo - E){1 - sin 6*+ sin6'_ cos0 - cos 61+ cos ) + w£'(cos6'+ - 1) + w(w - E){cos6^ - 1) + 



E is the positron energy and we have assumed that the polarization direction is along the x axis (see Fig. 
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FIG. 2. Spatial azimuthal distribution of a pair created by 100 MeV photon 

We also note that for exactly coplanar events (that is, for ^ = tt) it is more probable to find a pair in a plane 
perpendicular to the polarization vector that in a plane parallel to this vector, but this tendency changes rapidly as 
(j) decreases. Repeating this calculation for different values of energies of the incoming ray yields similar results. 

Since one of our goals is to use the azimuthal distribution in numerical simulations of the performance of a gamma 
ray detector, it is very convenient to parametrize the distribution surface in terms of smooth functions of (0, ■0) . 
This parametrization is done as follows. We first select coplanar events, cf) — tt. Instead of using a fitting function 
A*{1 — A*cos^i/;), used in previous works [0, we choose a form /(tt,?/;) — /^/2 sm^ i/) + /qCOS^V (note that the 
asymmetric ratio can be computed by the quotient ///7r/2)- As we move off this coplanar situation, /o and /7r/2 
become functions of 0, that is: 

/(</>, V') = L/2{^) sin" ^ + fo{^) cos2 . (6) 

Thus, to parametrize the entire surface we simply evaluate the azimuthal distribution for two values of ■0, namely, 
and 7r/2 . 

Since both fo{4>) and f^/2{4') ^re functions that rapidly vary when approaches tt, it was necessary to adjust the 
functions in two ranges of 0: (I) < < 3.05 rad. (II) 3.06 rad < < tt , whereas in the small range 3.05 < < 3.06 
we extrapolate the two fitting functions until the intersection point is reached (this minimizes the errors in the above 
mentioned range). 
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TABLE I. Parameters for the Lorentzian fitting /o (<?!>) and /,r/2(0) 



Lorcntzian 




/o(0) 








/7r/2(0) 






Energy 


yo 


A 






yo 


A 






100 MeV 


6.0 


1.03 


0.089 


3.118 


3.4 


1.84 


0.101 


3.146 


200 MeV 


8.16 


1.76 


0.0598 


3.130 


4.7 


2.7 


0.064 


3.148 


500 MeV 


10.8 


3.06 


0.036 


3.138 


6.3 


4.8 


0.034 


3.150 


700 MeV 


12 


3.6 


0.029 


3.139 


6.9 


6.0 


0.026 


3.15 


1 GeV 


13 


4.3 


0.0248 


3.141 


7.8 


7.8 


0.018 


3.15 



TABLE IL Fit for the parameter of fo{4>) function. 



Parameter 



Function 



yo 

A 



a\nE - b 
a\nE - h 
a + h/E + c/E^ 
a^h/E ^ c/E^ 



2.98 ±0.06 
1.41 ± 0.08 
0.015 ±0.001 
3.143 ±0.001 



7.7 ±0.4 
5.6 ±0.5 
9.5 ±0.6 
-2.7 ±0.2 



(-2.2 ±0.1)10'' 
(2 ± 1)10^ 



In region I we used Lorentzian functions of the form 



/(</-) = 2/0 



2Auj 



(7) 



(8) 



whereas for region II the best fitting function was found to adopt the form: 

/((/)) = a + dtan {h(f) ± c) . 
Table | and || show the results obtained. 

It is interesting to observe that for region I the difference between the calculated distribution (Fig. Ill) and the 
parametric form given by (||) with the parameters of Table || is less than 5%. Note also that the coefficients themselves 
vary smoothly over the energy range. It is thus possible to give a parametric form of these coefficients as functions 
of the incoming energy. In Tables |l| and [II we give the explicit forms of jo{E) , lj{E) , A{E) and Xc{E) associated 
with fo{4>) and f-K/2{4') with appropriate fitting functions. 

The parametric description of Region II is given below in Table It follows from this table that the parameters 
associated with the tangent function are approximately the same for /o((/>) and do not differ by much for /,r/2(</') ■ 
Thus, if we consider that this part of the azimuthal function is very small when compared with the range close to tt , 
we find it reasonable to take an average of these constants for the entire energy range. 



IV. DETERMINING GAMMA RAY POLARIZATION WITH THE ASYMMETRIC RATIO 

From the azimuthal probability distribution we define the asymmetric ratio (AR) as follows. We first select the 
azimuthal distribution at (j) — tt (coplanar events), two azimuthal angles "0 = 0, 7r/2 (parallel and perpendicular to the 
polarization direction, respectively) and numerically integrate the distribution around those values of ((/), V') in the 



TABLE TIL Fit for the parameter of /^/2(<?!>) function. 



Parameter Function 



yo alnE-b 1.85 ± 0.07 5.1 ± 0.4 

A alnE-b 1.3 ± 0.1 (6.6 ± 0.2)10"^ 

uj a + b/E + c/E^ 0.008 ±0.002 12.1 ± 0.9 (-2.8 ± 0.8)10* 

Xc 3.149 
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FIG. 3. Asymmetric ratio for several gamma - ray energies versus azimuthal resolution 

The reader should not conclude from the above paragraph that one should build a poor angular resolution detector 
(with great savings for the project) to measure polarized gamma rays. The analysis was done under the assumption 
that there was only one incoming polarization direction. If one has astronomical sources that produce two or more 
polarization directions one obviously needs higher resolution to resolve these directions. 
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FIG. 4. Comparison of three asymmetric ratio for 100 MeV incoming energy using different selection criteria. (Cross) without 
selecting energy; (Box) selecting particles with energy between 45 and 55 MeV; (Diamond) selecting between 40 and 60 MeV. 

The values of AR obtained in Fig. ^ can be improved if one selects particles with some specific values of energy. In 
particular, if one selects particles with energies close to equipartition of the incoming photon energy one can improve 
the value of the asymmetric ratio. 

In Fig. 1^ we compare three AR for incoming 100 MeV gamma rays. The lowest curve is obtained by selecting all 
possible energies, the second curve is obtained by selecting particles with energies w/2 ± O.OSw, whereas the last one 
has energies in the range w/e ± O.lw. 

Note that when one carefully selects the energy range of the electron - positron pairs, one improves the values of 
AR achieving a maximum of approximately 1.4 at 0.65 rad. 

On the other hand, this improvement of 10% is obtained at the expense of discarding more than 80% of the produced 
pairs. Thus, it would appear that to determine the AR it is more reasonable to simply integrate over the whole energy 
scale. 



V. SUMMARY AND CONCLUSIONS 

In this work we have explicitly computed the differential cross section for pair production and the azimuthal 
distribution of the produced pairs. We have also given a parametric form of this distribution in the different regions 
of interest. This parametric form is specially useful when implementing a Monte Carlo routine. Wc have then shown 
that the Asymmetric Ratio yields, in principle, a method to detect the polarization direction of the incoming photon. 
It is interesting to note that the azimuthal distribution surface changes rather abruptly for different values of (j) near 
coplanarity. The AR function reflects this behavior. It starts at a value of 0.8 for coplanar events and reaches a 
maximum value above 1.25 when the angular distribution is about 0.5 rads. The graph of AR explains and resolves 
the discrepancy between different authors who used this function to detect polarization direction. 

As we mentioned before the amplitude of the AR can be improved if one selects the energy range of the produced 
pairs. However, this improvement of approximately 10% is obtained at the expense of discarding more than 80% of 
the pairs. It then appears that unless the intensity of the incoming gamma rays is very high, it is more convenient to 
use the entire energy range of the outgoing pairs. 

The main conclusion one draws from the AR is that it is not necessary to have a detector with very good angular 
resolution to determine whether or not the incoming rays are linearly polarized. On the other hand, one needs very 
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good resolution to estimate the direction of the polarization vector. Thus, any 'reasonable' resolution that enables 
to determine the direction of the incoming photon Q should be adequate to then decide whether or not he photon is 
linearly polarized. The better angular resolution of the detector, the more accurate the determination of the incident 
direction and polarization vector of the gamma ray. 

Finally, a question we would like to address is whether or not it is possible to construct a detector capable of 
measuring the incident polarization of gamma rays produced by astrophysical sources. Preliminary results obtained 
with a monte carlo code simulation (where it was assumed that all events are coplanar) show that it is possible 
to determine the gamma ray polarization with an optimal geometrical setup . Future work will include the 
parametric form of the azimuthal distribution here discussed for a more realistic monte carlo simulation. We will also 
analyze alternative geometries of the measuring apparatus to minimize the errors introduced by multiple scatterings 
of the produced pairs. 
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